Evolution of spectral properties along the 0(6)-U(5) transition 
in the interacting boson model. I. Level dynamics 
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We investigate the evolution of quantal spectra and the corresponding wave functions along the 
[0(6)— U(5)]DO(5) transition of the interacting boson model. The model is integrable in this regime 
and its ground state passes through a second-order structural phase transition. We show that the 
whole spectrum as a function of the Hamiltonian control parameter, as well as structures of all 
excited states, exhibit rather organized and correlated behaviors, that provide deeper insight into 
the nature of this transitional path. 
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I. INTRODUCTION 

Properties of the interacting boson model (IBM) [IJ in 
transitional regimes between various dynamical symme- 
tries have been extensively studied mainly in connection 
with zero-temperature quantum phase transitions 0, HI • 
In any of such transitions, the structures of the ground 
state and few low-lying states change abruptly (for the 
system size tending to infinity) at a certain critical point, 
located on the way between the two dynamical-symmetry 

limits, see, e.g., Refs. HHUHBH- This behavior finds 
experimental evidence (in a finite- N approximation) in 
observed variations of nuclear shapes in some isotopic or 
isotonic chains of nuclei. 

The IBM phase transitions are of the first order, except 
the isolated point of a second-order transition, which is 
located at the intersection of borders between spherical 
and deformed, and between prolate and oblate shapes 
in the parameter space 0,H|,Ei- To pass this point, one 
commonly starts from the 0(6) dynamical symmetry and 
proceeds to U(5) via the line of unbroken 0(5) dynam- 
ical symmetry. The deformcd-to-spherical second-order 
phase transition on this path manifests itself as a non- 
analytic but continuous change of the ground-state de- 
formation, in contrast to the discontinuous changes ob- 
served along the other (even infinitely close) transitional 
paths. This type of phase structure of the parameter 
space agrees with the classical Landau theory of thermo- 
dynamic phase transitions, that is applicable at zero tem- 
perature if the role of thermodynamic variables is taken 
by the model control parameters 0, II , and with catas- 
trophe theory 0, HJ . 

The above-mentioned [0(6)— U(5)]dO(5) transitional 
path differs from the others also in that it does not de- 
stroy the integrability of the Hamiltonian. Due to the 
underlying 0(5) dynamical symmetry [To| . the integrals 
of motions along the whole transition form a complete set 
of commuting operators and the Hamiltonian eigenprob- 
lem can be solved analytically [llj, [l2| . This was used 
for an explicit calculation of some second-order phase- 
transitional observables Recent studies of the 0(6)- 



U(5) transitional path were also based on the concepts of 
the E(5) critical-point dynamical symmetry [lj, [l5j and 
the quasidynamical symmetry [l6| . 

The ultimate mechanism that is on the deepest level 
responsible for the occurrence of ground-state phase tran- 
sitions of various orders in quantum many-body systems 
remains unclear. The distinction between the IBM first- 
and second-order phase transitions, e.g., was shown [t| 
to be connected with different densities of unavoided en- 
ergy crossings (branch points) in the complex-extended 
parameter space, which in the N — > oo limit accumu- 
late infinitely close to critical points on the real axes (in 
analogy with similar behaviors of complex zeros of parti- 
tion functions in thermodynamic phase-transitional sys- 
tems). However, many questions — among them the role 
of integrability in the process of dynamical-symmetry 
breaking — still remain open. 

The present work contributes to the mapping of this 
relatively new territory of physics by studying in de- 
tail various spectral observables associated with the in- 
tegrable phase-transitional path in the IBM. In partic- 
ular, we investigate the evolution of energies and wave 
functions of individual Hamiltonian eigenstates with zero 
angular momentum along the whole [0(6)— U(5)]dO(5) 
line. It is shown that this transitional class of IBM ex- 
hibits rather peculiar features. 

We will combine two totally different, but mutually re- 
lated general approaches: (i) the theory of level dynam- 
ics, initiated by Pechukas and Yukawa [l7|. also known as 
the dynamical Coulomb-gas analogy, and (ii) the semi- 
classical theory of quantal spectra, represented by the 
Gutzwiller and Berry- Tabor trace formulas [l^. Results 
obtained by applying both these approaches will be pre- 
sented in two parts: approach (i) is discussed in the 
present article (Part I) , which gives numerical results on 
level dynamics, while approach (ii) will be used in the 
following article HH (Part II). 

The Pechukas- Yukawa theory describes the dynamics 
of individual levels and the interaction matrix elements 
via a set of coupled differential equations, where the vary- 
ing Hamiltonian control parameter plays the role of time. 
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It enables one to understand the evolution of spectral ob- 
servables with the control parameter (including eventual 
phase transitions) in a more intuitive way, using the par- 
allel with a classical ensemble of charged particles moving 
in one dimension. 

The trace formulas, on the other hand, describe a snap- 
shot of the energy spectrum at each fixed value of the 
control parameter ("time") by expressing the quantum 
density of states (as a function of energy) through prop- 
erties of periodic orbits in the classical limit of the sys- 
tem. Since both methods (i) and (ii) translate the orig- 
inal problem to the classical language, they often pro- 
vide deeper understanding of specific behaviors observed 
on the quantum level. Also in our case, the most sig- 
nificant features of the IBM in the [0(6)-U(5)]dO(5) 
transitional regime will be elucidated by both kinds of 
classical concepts involved in the above approaches. 

The plan for this part of the paper is the following: In 
Section [H] we will briefly describe the quantum Hamil- 
tonian under study, its integrals of motions and phase- 
transitional features. Section II I II presents numerical re- 
sults on the level dynamics and their interpretation in the 
framework of the Pechukas- Yukawa theory. The accom- 
panying changes in the structure of wave functions are 
then discussed in Section llVl Finally, Section Ivl contains 
partial conclusions of this part of the paper. 



II. QUANTUM HAMILTONIAN 

The interacting boson model pj describes shapes and 
collective motions of atomic nuclei in terms of an ensem- 
ble of N interacting s and d bosons with angular mo- 
menta and 2, respectively. To analyze the evolution of 
properties of this model along the 0(6)-U(5) transitional 
path, we adopt the Hamiltonian 



H( V ) = a 



(1) 



where the dimensionless control parameter 77 € [0, 1] 
changes the proportion of both competing terms and 
drives the system between the 0(6) (77 = 0) and U(5) 
(77 = 1) dynamical symmetries. The operator fid — 
[d) ■ d) = N — n s represents the (i-boson number, while 
Q = Q { 2) = [std+dtg](2) stands for the 0(6)-U(5) 
quadrupole operator. The energy scale is set by an arbi- 
trary factor a, which in the following will be fixed at the 
value a = 1 MeV. 

Hamiltonian (JIJ is a special case of a more general 
Hamiltonian of the same form, but with the quadrupole 
operator given by = Qo + xl^d}^, where x € 

[— ^ is an additional control parameter. Eq. Q), 
where % = 0, can be decomposed [2fJ into a linear combi- 
nation of Casimir invariants corresponding to the 0(6), 
0(5), 0(3), and U(5) algebras, with no admixture of 
SU(3), SU(3), and 0(6) invariants, i.e., describes the 



[0(6)— U(5)]dO(5) transitional line in the extended Cas- 
ten triangle |6J. 

The above Hamiltonian can also be rewritten as 

H( V ) = (1 - r,)H(0) + V H(1) = H + VV , (2) 

which is the form well known from various studies of 
quantum phase transitions. Assuming a = 1, we obtain 



(3) 
(4) 



The evolution of Hamiltonian @ with 77 can be treated in 
a perturbative way since H(rj + 577) = H (7/) + 5rj V . Note 
that the powers of N in denominators of Eqs. Q), (0 
and Q guarantee convenient scaling of the Hamiltonian 
with variable boson number N 1. 

It can be easily shown that for Hamiltonian J5J the 

ground-state average (V) v = (i>i(r))\V\ipi(r))) = 
[where Ei(jf) and |V>i(t7)) are the ground-state energy 
and wave function, respectively] is a nonincreasing func- 
tion of 77. Therefore, if V is nonnegative — as in our spe- 
cific case, see Eq. (0J — then an instantenous satisfaction 
of (V) Vc = at some critical point ?7 C implies that the 
average gets fixed for all 77 > ?7 C , freezing both the en- 
ergy and wave function of the ground state. At this 
point, the system may exhibit (for N — * 00) a ground- 
state phase transition of order k > 2. If the second 
derivative of energy changes discontinuously from a value 
d ^2 V) = < at 77 = ?7 C _ to zero at 77 = ?7 C+ , the 

transition is of the second order. Higher-order transi- 
tions 0, El would require additional constraints, namely 

=0forfc<K. 

ar} K ir/c- 
It is not difficult to see that for the specific Hamil- 
tonian in Eq. the ground-state average of V indeed 
interpolates between a positive value at 77 = and zero 
at 77 = 1. However, the phase-transitional scenario is 
generically allowed only in the limit of infinite Hilbert- 
space dimensions, thus N — > 00, when the ground- 
state energy as a function of 77 may acquire nonanalytic 
character. The asymptotic critical point is located at 
77c = § = 0.8. At this point, the deformed ground- 
state configuration, given by a mixed-boson condensate 
\ipi) cx (s^ + /3 gs d,Q) N \Q), changes into the pure s-boson 
condensate, \tpi) oc (sT^O), characterizing the spheri- 
cal U(5) phase. In the left vicinity of the critical point 
the ground-state "deformation parameter" [3 gs drops to 
zero as /3 gs oc ^/rf c — 77 7] and the corresponding value 
of (V)rj behaves according to N — > 00 asymptotic for- 
mula (y)n oc (77c — 77) 9]. Thus both f3 gs and {V}^ can be 
considered as order parameters describing a second-order 
quantum phase transition, k = 2, with critical exponents 
\ and 1, respectively. 

The limits H(Q) and H(l) of Eq. JTJ posses the 0(6) 
and U(5) dynamical symmetries, respectively. Since the 
dynamical-symmetry Hamiltonians are constructed using 
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solely observables "in involution" (the Casimir invariants 
of the respective algebraic chain), they are always inte- 
grable .21] . Moreover, because the 0(5) dynamical sym- 
metry underlying both 0(6) and U(5) limits is not broken 
in the transitional regime, the integrability of Hamito- 
nian JJJ is preserved for all values of rj [2H 122^ . Indeed, 
one can find five mutually commuting integrals of mo- 
tion, the same number as the dimension of the classi- 
cal configuration space (given by two geometric parame- 
ters and three Euler angles 23]). Four of these integrals 
can be associated with the following quantum numbers: 
energy Ei(rj) given by the Hamiltonian H(t)), squared 
angular momentum 1(1 + 1) represented by L 2 (where 
L = [Sdy- 1 '), its projection m determined from L z , and 
the seniority v defined through the v(v + 3) eigenvalue of 
the 0(5) Casimir invariant 



C 2 [0(5)] = -(L-L) + 2(T 3 -T 3 ) 
5 



(5) 



(where T 3 = [cftd]w). The fifth integral of motion, 
connected with the so-called missing label h& of the 
0(5) DO (3) reduction, is not given explicitly, but its ex- 
istence is guaranteed by the fact that there must be five 
independent commuting operators in the complete set, 
so the Hamiltonian (which is made of four of them) com- 
mutes with the fifth one |2l|. Note that in this paper 
we will only consider the set of states with zero angular 
momentum, / = 0. 



III. EIGENVALUE DYNAMICS 
A. Pechukas- Yukawa equations 

Drawing the dependence of all individual level ener- 
gies Ei(rj) for Hamiltonian © on the control parame- 
ter, one obtains a picture containing n continuous curves 
that resemble trajectories Xi(t) of an ensemble of parti- 
cles in one dimension. The motion of levels is described 
by a set of Hamilton-type first-order differential equa- 
tions associated with a gas of particles interacting via 
two-dimensional Coulomb force: 



cPEj 
dr] 2 



= 2£ 



\Vy\ 2 

Ei — E.j 



(6) 



(analogous to 



dt 2 2ire ^-'j(¥ : i) Xi 

to the ordinary gas dynamics, however, the 
charge" \Vij\ 2 = \(ipi(r])\V\^j{r)))\ 2 <-> = Q i3 can- 
not be factorized and varies as the "time" 77 «-> t elapses. 
Thus the product charges (alias interaction matrix ele- 
ments) are also dynamical variables, subject to specific 
evolution, and the system's phase space is larger than 2n. 
Besides Eq. © we have 
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In contrast 
product 



dVjj 
dr] 



Vu - V, 
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Ei + Ej — 2Ek 
(Ei — Ek)(Ej — Ek) 



V ik V ki 
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FIG. 1: Spectrum of Hamiltonian (1) with N = 40 as a func- 
tion of r] for 1 = levels with all seniorities. The vertical and 
horizontal rectangles are expanded in the upper right and 
lower panels, respectively. The seniority is assigned to several 
levels in the lower panel. 



for i 7^ j, and 



dEj 
dr] 



(8) 



(7) 



Eqs. iJBJ 10) are equivalent to the well-known 
Pechukas- Yukawa set of equations ^tJj although we use 
here a slightly different form than the one usually found 
in textbooks |2J|. The system described by these equa- 
tion is deterministic and even integrable. If all energies 
and interaction matrix elements are known at a single 
point rj (for instance 77 = 0), the equations determine 
Ei's and Vij's for all other 77 values. 

Since the product charge |Fy| 2 in Eq. JBJ) is non- 
negative, the levels never touch each other unless their 
mutual interaction completely vanishes. In absence of 
symmetry-dictated zeros of the interaction matrix, the 
coincidence of simultaneous convergences \V i:j \ 2 -> and 
Ei + \ — Ei — > is extremely unlikely, which gives rise to 
the well-known "no-crossing" rule for level energies. The 
presence of symmetries, however, induces the disappear- 
ance of Vij's for certain sets of states which, therefore, 
can cross. In case of Hamiltonian (JTJ, this concerns levels 
with different values of angular momentum I and levels 
with different seniority v. 



B. Level bunching around E « 

In Figure ^ we show the dynamics of all levels with 
I = along the 77 E [0, 1] path between the 0(6) and 
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U(5) dynamical symmetries of Hamiltonian Q. The cal- 
culation was performed by numerical diagonalization of 
the Hamiltonian for N = 40 bosons. One can observe nu- 
merous level crossings, particularly in the region around 
E rs (the horizontal rectangle, expanded in the lower 
panel), which is a consequence of the unbroken 0(5) dy- 
namical symmetry of the system. Indeed, the seniority 
quantum numbers, as marked for few levels on the right- 
most side of the lower panel, differ for any pair of levels 
that cross at some point. We will see below that the 
crossings disappear after separation of levels with differ- 
ent seniorities into several figures. 

The pattern of consecutive compressions and dilutions 
of the spectrum in the region around E ss is one of 
the most apparent attributes of Fig. ^ ( see the lower 
panel). A striking feature of this pattern is the regular 
sequence that characterizes the total numbers of levels 
involved in individual bunches: when descending from 
r\ = 0.8 to w 0.4, the sequence goes like 1, 2,3,4, ... At 
first, the different seniority states seem to cross exactly at 
the same point (within the available numerical precision) , 
but with rj descending below 0.65 the higher seniorities 
get increasingly out of focus, and the bunching pattern 
becomes more and more diffuse. Nevertheless, the struc- 
ture of alternating clusters and gaps extends over a wide 
range rj S [0.3,0.8]. Secondary "interference" patterns 
are also visible at other energies (see the vertical rect- 
angle of Fig. n extended in the upper right-hand-side 
panel), but these are much weaker than the main one. 

The energy E f=a 0, where the bunching pattern ap- 
pears, is significant because it corresponds to the local 
maximum at f3 = of the classical potential P, |2^| cor- 
responding to Hamiltonian The bunching of levels 
thus develops just at the value of energy where the clas- 
sically accessible range of the deformation parameters, 
/3 £ [Anin> Pmax]t extends due to /3 min becoming zero. 
The connection of the bunching pattern with the IBM 
classical dynamics will be elaborated in Part II of this 
contribution [Igj . 



C. Shock- wave scenario 

Figure |21 demonstrates that the level bunching pattern 
can be deconvoluted by separating states with different 
seniorities. Here we show the level dynamics for v = 
(this set includes the ground state) and v = 18, with 
the boson number N = 80. Clearly, the v — levels 
in panel (a) form a smooth flow with a "shock wave" 
propagating from the top of the spectrum (at 77 = 0) 
to the ground state (at rj = 0.8). The mutual distances 
AEi = Ei + i — Ei of individual v = levels as functions of 
r\ are shown in Figure Ola), where we can clearly identify 
points of the closest approach of neighboring states as 
the wave propagates through the ensemble. Note that 
due to the energy denominator in Eq. (jSJ, a minimal 
spacing of levels tends to induce maximal "force" acting 
on the relevant levels, which is basically the mechanism 
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FIG. 2: Spectrum of Hamiltonian (1) with N = 80 for I = 
levels with seniority v — (panel a) and v = 18 (panel b). 
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FIG. 3: Distances of neighboring levels from Fig. |21 
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FIG. 4: Branch points of Hamiltonian Q with AT = 20 for 
Z = states with seniorities v = (panel a) and w = 6 (panel 
b). 



the corresponding levels, which does not typically hap- 
pen (for levels with different symmetry quantum num- 
bers). On the other hand, if a given branch point is not 
on, but sufficiently close to the real axis, one observes 
an avoided crossing of the relevant levels at the corre- 
sponding value of rj. Remind that a sequence of such 
avoided crossings is significant for the "shock wave" . A 
cumulation of branch points in infinitesimal vicinity (for 
N — > oo) of the critical point -q c was recently shown 
to constitute the essential triggering mechanism for the 
fBM quantum phase transitions of both orders. 



For each seniority, there are altogether 



n(n— 1) 



complex 



conjugate pairs of branch points, where n is the dimen- 
sion of the given seniority subspace. Because of numer- 
ical constraints, we can only show results for moderate 
dimensions that correspond to the lower boson number 
TV = 20. As can be seen in Fig.QJ branch points for both 
seniorities form rather regular patterns, fn the v = 
case (panel a) we notice a chain of points at rj < 0.8 that 
approach close to the real axis. These points clearly cor- 
respond to the sequence of avoided crossings shown (for 
a higher boson number) in Fig. HJa) . With increasing se- 
niority, the pattern gets more and more separated from 
the real axis (see the example in panel b), which results 
in a weakening of level interactions, as observed (for dif- 
ferent values of N and v) in Fig. E{b). Note that such 
an organized behavior of branch points is a privilege of 
only the [0(6)— U(5)]dO(5) transitional class, where the 
separation of seniorities is possible (cf. Ref. 



that keeps the wave moving. This is also why the wave 
initiates in the upper (densest) part of the spectrum (at 
rj = 0, the distance of nearest levels linearly decreases 
with i, while at 77 = I it is constant, cf. Fig- • 

On the other hand, the dynamics of the v = 18 levels, 
shown in panels (b) of both figures, exhibits much weaker 
interactions. The flow in Fig. E^b) looks almost laminar 
and the minimal distances in Fig. [2b) (still disclosing in- 
teractions) are about twice larger than in the v = case. 
It can be checked that the weakening of level interactions 
proceeds gradually as v increases. 

The shock- wave interpretation of Fig. HJa) is particu- 
larly appealing if used as a tentative reasoning for the 
ground-state phase transition at ry c = |. It seems that 
this transition results from a highly ordered sequence of 
structural changes that propagate from upper to lower 
parts of the spectrum and terminate at the ground state 
just at the critical point. This mechanism, however, 
needs to be verified by an analysis of wave functions and 
will be further discussed in Sec. IIVI 

Closely related to the regular evolution of level ener- 
gies is the organized pattern of the Hamiltonian branch 
points in the complex plane of parameter rj 25]. It is 
shown for N = 20 in Figure 01 for (a) v = and (b) 
v = 6. Branch points are places in the complex-extended 
parameter space where two (or more) Hamiltonian com- 
plex eigenvalues become degenerate j2g ■ A branch point 
located on the real 77-axis would imply a real crossing of 



D. Focal point and spectral invariant 

A more detailed view of Fig. |2fa) discloses that almost 
all v = levels on the 77 = side (except perhaps few 
ones on the top of the spectrum) point to a virtually 
sharp focus on the r) = 1 side. Indeed, an unperturbed 
evolution (with no mutual interactions between levels) 
would lead to crossing of individual lines at the point 
(t],E) — (1,2)) which we call an (approximate) focal 
point of the [0(6)-U(5)]dO(5) transition. 

From Eq. JHJ) we see that (r/f , Ef) will be a focal point 
of Hamiltonian 10) if (ipi(0)\H(Tjt)\ipi(0)) = E{, so in our 
particular case we have 

(^(0)1^(0)) « ^ , (9) 

where \4>i{0)} are the Hamiltonian eigenvectors with v — 
at rj = 0. This means that the average number of d- 
bosons in individual 0(6) eigenstates with zero seniority 
stays nearly constant across the whole spectrum. Fig- 
ure |Sf a), where the average is shown explicitly along 
the whole 77 € [0,1] path for all v = levels with N = 80, 
supports this rule; see the rj = limit (graphically it 
is difficult to distinguish, whether the convergence of all 
curves to the y point is exact or not, but numerical val- 
ues indicate that it is only approximate). With a lower 
precision, the validity of the above "spectral invariant" 
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FIG. 5: The average number of d-bosons for v — (panel a) 
and v = 18 (panel b) states with I — and N = 80. 

can be extended to higher seniorities, but with increasing 
v there are more and more upper states that do not fit, 
see Fig.|5fb) that shows (rid)i for the v = 18 levels. 

Even in the SU(3)-U(5) and SU(3)-U(5) transitions, 
when Q in Hamiltonian is replaced by Q x with 

X = ~ -y or x = +y and the seniority is not con- 
served, one finds a similar approximate invariant, namely 
(^(0)1^1^(0)) re M;, where |^(0)) represent the SU(3) 
or SU(3) eigenvectors (the complete I — spectrum for 
these transitions can be found in Ref. §]). 

Note that we first detected these invariants geometri- 
cally, from the focal points. The impact of such invariants 
on the level dynamics can be enormous since focal points 
represent an essential condition for the initial compres- 
sion of the Coulomb gas, which results in stronger in- 
teractions between levels. This compression triggers the 
formation of the "shock wave" in the densest part of the 
spectrum, see Sec. MI CI Therefore, the existence of an 
initial (exact or approximate) focal point may belong to 
the main causes that eventually lead to a phase transition 
at some point. 



E. Finite-iV phase transitions 

Although we are dealing here with the spherical- 
deformed transition induced by varying parameter 
r\ in Hamiltonian (JTJ, one should realize that the 
[0(6)— U(5)]dO(5) transitional path itself coincides with 



the separatrix between prolate and oblate deformed 
phases 0, EJ • The prolate-oblate first-order phase tran- 
sition for N — » oo at any fixed value of r\ £ [0, |) can 
be induced by varying parameter x i n the generalized 
Hamiltonian of the form with Q replaced by Q x . 

It is argued in Ref. ^3] that the 0(6) dynamical sym- 
metry represents a special point of the IBM phase dia- 
gram where a discontinuous prolate-oblate change of the 
ground-state structure can be observed even for finite bo- 
son numbers. Indeed, if one explicitly includes the 0(5) 
Casimir invariant JSJ into the Hamiltonian with a coeffi- 
cient such that the v = ground-state at (ij, x) = (0, 0) 
becomes degenerate with the lowest states of other senior- 
ities, a crossing of the ground-state configurations will 
occur for any value of N when passing the 0(6) point in 
the ^-direction. 

It is easy to show that this mechanism can be extended 
to the whole r\ £ [0, |) transitional region. The basic 
trick — the fact that levels with different seniorities can 
be made degenerate — remains the same. After subtract- 
ing the component corresponding to the 0(5) Casimir 
invariant | 2f| from the general x-dependent Hamiltonian 
of the form |JT|I. one arrives at the expression 

H'(V,X) c< ^ |(Q X • Q x ) - ^2 [0(5)]} + ^n d , 

(10) 

that exhibits the desired property: For any fixed value of 
■q and any finite boson number N, the ground state as a 
function of x changes discontinuously at x.c — 0, which 
thus defines the critical point for a finite- iV prolate-oblate 
phase transition. 

We therefore extend the region of possible finite-A r 
phase transitions in the IBM phase diagram to the whole 
prolate-oblate separatrix. The key feature needed for 
this generalization of Ref. is the integrability of the 
[0(6)— U(5)]dO(5) Hamiltonians. Note, however, that 
phase transitions at finite dimensions, induced by un- 
avoided crossings of levels involving the ground-state, 
represent a rather nongeneric kind of behavior, which 
moreover is not robust enough to survive at finite tem- 
peratures. Indeed, if the temperature increases from zero 
to an infinitesimally small value, nonzero populations of 
both levels result in a smooth dependence of the free en- 
ergy on the control parameter, and the phase-transitional 
behavior is washed out. 



F. Bulk properties of the spectrum 

It is clear that the strongest influence on a given 
level comes typically from its neighbors at the places of 
avoided crossings. Besides these binary interactions (in- 
volved in the shock- wave propagation), there exists also 
a component of the total force acting on each individual 
level that originates from the bulk of the whole ensemble. 
In this subsection, we will consider two global measures 
of this bulk component. 
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FIG. 6: The dispersion I1H of the spectrum and the kinetic 
energy from Eq. 1131 for v = and 18 states, corresponding 
to I = and N = 80. 



First, we consider the overall compression of all levels, 
represented by the energy dispersion (squared "spread" ) 
of the spectrum, A|, = rjEiC^i ~ ^) 2 > where E = 
— £"i is a center-of-mass energy. A straightforward 
calculation yields the expression 



A 2 _ 
L± E — 



Tr 2 H 



2, 7 



Tr(H V) TiHqTiV 



TtV 2 



Tr 2 V 



(11) 



which shows that the spectral dispersion is a quadratic 
function with a minimum at 



Vo = 



nTijHoV) - TrgpTrV 
nTrV 2 - Tr 2 V 



(12) 



For 77 ps 770, the strengths of both terms i?o and r/V of 
Hamiltonian J5J are comparable, so that the strongest 
effects of mixing take place in the surrounding region. 
For 77 3> 770 or 77 <C 770, on the other hand, the spectrum 
just blows up, the Hamiltonian being dominated by r/V. 

For (a) v = and (b) v = 18 subsets of the spectrum 
with N = 80, the function (|11|) is shown in the two left- 
most panels of Figure We see that the v = levels are 
maximally compressed at 770 ~ 0.56, i.e., in the region 
just before the phase transition. For higher seniorities, 
the minimum moves towards r] c = |. Let us note that a 
similar conclusion can be made for x ^ 0, when of course 
the seniority is not conserved and the contribution of all 
/ = levels must be summed up. For )( = ±y, for 
instance, the energy dispersion forms a sharp minimum 
directly at 770 « 0.8. 

The second quantity we will use here to characterize 
the bulk component of the force is the total product 



charge Q = Y,i>jQi] 



E l>7 l^l 2 - 



the sum 
2 ^ 



(13) 



T 



It is related to 



which is an integral of motions of the Pechukas- Yukawa 
model, known as the total energy [24|]. Since -^E = 0, the 
second term that represents the potential energy V = Q 
is at any value of 77 just a complement of the first, kinetic 
term T. For r\ 3> ?/o and 77 <C 770 [assuming for a while 77 S 
(— 00, +00)], the eigenbasis of -ff(?7) virtually coincides 
with the eigenbasis of V so that TrV" 2 rs Y, v i = 27 " 
and V ~ 0. In these regions, the gas just freely expands. 
On the other hand, around 770 the kinetic and potential 
terms in Eq. (|13l) are comparable and the interaction may 
generate nontrivial effects. 

The kinetic energy from Eq. i|13|) for levels with v — 
and v = 18, respectively, is shown in panels (c) and (d) 
of Fig. [f)| In both cases, we observe a minimum of T 
very close to the critical point; for v — the minimum is 
located at 77 1=3 0.67. This means that V = Q is maximal 
at the same place, implying the strongest overall strength 
of level interactions. For higher seniorities, the minimum 
gets shallower and moves towards ?7 C . 

We saw that both the compression of the spectrum 
and total interaction strength are maximal in the region 
of control parameters around 770 which immediately pre- 
cedes the phase transition at ?7 C . Conversely, Eq. ((T^l) 
yields a reasonable rough estimate of the parameter range 
of a general Hamiltonian J5J) where eventual phase tran- 
sitions may be located. 



IV. EIGENSTATE DYNAMICS 

Besides dynamics of individual Hamiltonian eigenval- 
ues Ei(rj), one can also analyze structural changes of 
the corresponding eigenstates 1 7/^(77)). These two as- 
pects of spectral evolution are mutually correlated, since 
the matrix elements Vij , that carry information on wave 
functions, belong to dynamical variables involved in 
Pechukas- Yukawa equations 

In Fig. [SI we have already seen the evolution of the av- 
erage number of <i-bosons, (nd)i, in the v = and v — 18 
eigenstates. This information is now supplemented by 
Figure where the 77-dependence of the whole distribu- 
tion Pi{rid) of rid is shown for selected N — 80 Hamilto- 
nian eigenstates, namely the v = states with i — 1, 10, 
20, and 30 (ordered with increasing energy). For I and v 
fixed, the probability 

Pi(n d )\ = ^2 K^^n-A^HV^))! 2 (14) 



for each 77 is determined as the projection of the state 
|t/^(t7)) onto the subspace of the U(5) eigenstates with 
rid equal to the given number. In the U(5) limit, the 
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FIG. 7: The distribution of the d-boson number rid in four 
I = v = eigenstates of Hamiltonian Q with N = 80 as a 
function of rj. The insets show the corresponding U(5) wave- 
function entropy. 



distribution is concentrated on a single value rid — 2i — 2 
(with zero seniority, the value of rid must be even), but it 
quickly spreads over a broad range of rid's as rj decreases 
from 1 to 0. 

Figure [7| shows four qualitatively different types of 
ways how this derealization proceeds: For the ground 
state, i = 1, the value of rid remains zero as far as 
rj > rjc, and then it suddenly increases (with decreas- 
ing rj), forming a ridge around the average that goes ap- 
proximately as {rid) i « \Jf)c — T), in agreement with the 
phase-transitional predictions; cf. Fig.|3{a). For excited 
states, the gradual spread of wave functions in rid can be 
compared to the propagation of waves on a string. The 
string is initially (at rj = 1 ; the "time" is now thought to 
go backwards) subject to an instantenous point pertur- 
bation and the resulting waves propagate in both rid = 
and rid = N directions asymmetrically. The pattern of 
wave propagations changes with i: for instance, the speed 
of the upper wave is lower for higher excited states. When 
the lower front of the wave reaches the rid — limit, it ei- 
ther gets reflected (this happens for lower excited states, 
see the i = 10 example) or stops there (for higher excited 
states, see i = 20 and 30 cases). 

It is interesting that the value of rj where the wave 
reaches the lower endpoint rid = coincides with the 
range where the "shock wave" affects the given level, see 
Figs. GJa) and[3fa). This can be checked for a larger 
set of levels in Fig.|5{a). The dependences of individual 
rid averages exhibit well-pronounced minima that corre- 
spond to the stopping or reflection of the lower wave front 
at rid = 0, and reasonably coincide with the moments of 
passage of the shock-wave. 



Also shown in the insets of Fig. [7| is the U(5) wave- 
function entropy, 

N 

^ U(5) =- J2 Pi(nd)lnPi(n d ) , (15) 

n d =0 

which measures the overall spread of the instantenous 
eigenvector \ipi(rj)) in the U(5) basis [2(|. Assuming a 
quasiuniform distribution of the ith state over a certain 
set of the fid eigenstates, one finds that the effective num- 
ber of components is given by n° e = expS^ 5 \ This 
number is approximately equal to half of the width (at 
a given value of rj) of the rid distribution corresponding 
to the respective level (taking into account that odd rid 
values are not populated for v — 0). 

As can be seen in Fig. [7| the widths of the rid distri- 
butions and the corresponding U(5) entropies grow with 
decreasing rj as far as the distribution touches the rid = 
limit (the level gets into the shock- wave region). Af- 
ter this point, the width and entropy stay approximately 
constant. If proceeding from the 0(6) side, i.e., return- 
ing to the forward direction of "time" , we can say that 
the process of localization of level i in the U(5) basis 
starts approximately when the shock wave hits the level. 
This supports and further specifies the shock-wave sce- 
nario described in Sec. IIII CI We must stress, however, 
that for excited states the transition to the U(5) struc- 
ture after passing the shock wave is only gradual. A 
sudden phase-transitional type of change is reserved for 
the ground state only. 

As indicated by the i > 1 examples in Fig. [7| the de- 
crease of the U(5) wave-function entropy exhibits some 
undulations, connected with quantum interferences of the 
amplitudes corresponding to populations of individual 
rid's. It is surprising that vertical coordinates of the main 
oscillations are about constant for the whole ensemble of 
states. This is demonstrated in Figure 00 where we show 
the U(5) wave-function entropy for all (a) v = and (b) 
v = 18 states (N = 80). Clearly, if one proceeds from 
state to state, the undulations are shifted in rj, but re- 
main at about the same levels of entropy. The result is 
a peculiar pattern of "plateaus" present in both panels 
of Fig. |S] (Let us stress, however, that these plateaus 
are only a visual effect appearing when all entropies are 
drawn in the same figure.) This hints at strong correla- 
tions in the structural changes of individual eigenstates 
after the passage through the shock-wave region. 

The most distinguished steps of the patterns in Fig. 
are the same for both seniorities. They correspond to 
the effective numbers of wave- function components equal 
approximately to n° ss 4.5, 7.5, 10, and 12. Note that 
the average derealization of a given state in a randomly 
chosen basis is for sufficiently high dimensions n given 
by n GOE ~ °- 48rl H3, which for the v = and v = 18 
subspaces yields typical saturation values of the wave- 
function entropy equal to Sqoe ~ 3 and ~ 2.7, respec- 
tively. (The largest U(5) entropies in the rj = limit 
slightly exceed the GOE values, but the latter provide 
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FIG. 8: The U(5) wave function entropy for all v — (panel 
a) and v = 18 (panel b) eigenstates of Hamiltonian Q with 
I = and N = 80. 

reasonable estimates of averages if all states are taken 
into account.) We see that the system of plateaus in 
Fig. |H1 disappear in noisy oscillations just below the re- 
spective GOE entropy values. 

Let us stress that no step-like structures are observed 
in cumulative plots of the U(5) wave-function entropy 
of all Z = states for the SU(3)-U(5) and SU(3)-U(5) 
transitions. The present correlated behavior is therefore 
connected solely with the integrable x = region. 

V. CONCLUSIONS 

We have studied dynamics of the I = energy 
levels and the corresponding eigenstates along the 
[0(6)— U(5)]dO(5) transition of the interacting boson 
model. Results of our numerical calculations were dis- 
cussed in the framework of the Pechukas- Yukawa model, 
which describes the evolution of quantal spectra as one- 
dimensional motions of an ensemble of classical particles. 
Treated in this way, spectral attributes for all values of 



the control parameter — including possible phase transi- 
tions at some critical points — result just from a specific 
"initial condition" , i.e., the set of energies and interaction 
matrix elements at a single arbitrary point r\. Of course, 
particularly tempting is to consider the whole spectral 
evolution along r\ G [0, 1] (and beyond) being predeter- 
mined by properties of the system in either of the two 
limiting dynamical symmetries. 

We disclosed cooperative and highly coherent behav- 
iors of the individual spectral constituents, i.e., level ener- 
gies and wave functions corresponding to various senior- 
ities. This may be generally linked to the integrability of 
the model in the present regime, namely to the possibility 
to separate seniorities, but we have to admit that some 
of the findings remain just plain observations. Further 
studies may still shed more light on how this all "comes 
about" . 

The most significant cooperative effect seems to rely 
on the "shock-wave" mechanism, that consists in an or- 
dered sequence of avoided crossing of levels in the re- 
gion around E w and the accompanying changes of 
eigenstates (Sees. lillCl and ITvjl . Triggered by an initial 
compression of the spectrum, the shock wave initiates 
in its densest upper part and propagates downwards to 
the ground state. The passage of the wave through a 
given state starts the gradual transfiguration of the state 
structure into the U(5) form. This mechanism provides 
a deeper insight into the process that eventually leads to 
the ground-state phase transition of second order. 

Among the other findings we highlight the following: 
(a) approximate focal points of IBM spectra in transi- 
tions to the U(5) dynamical symmetry (Sec. IIIID|) . (b) 
possible finite- ./V prolate-oblate phase transitions along 
the whole [0(6)-U(5)]dO(5) separatrix (Seo. lTTTEll . (c) 
extremes of spectral "bulk observables" in the region im- 
mediately preceding the phase transition (Sec. lIIIF|) . (d) 
highly correlated changes of consecutive eigenstates lead- 
ing to "plateaus" in the cumulative plot of U(5) wave- 
function entropies (Sec. lIV|) . 

In the following part of our article we wm f°~ 

cus on the interpretation of the E « pattern of level 
bunchings (Sec. IIII B)l within the semiclassical theory of 
quantal spectra. 
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